
**P-values for Table 4**

use "$data/merge_bl_fu", clear

*Initial data steps
	keep inc_asset_aggregate3_fu co_outcomes_aggregate_fu outlook_life_aggregate_fu ///
		hh_decision_pca_fu inc_asset_aggregate3_bl outlook_life_aggregate_bl ///
		hh_decision_pca_bl female d_lleligible bl_bus_owner fu_dummy* EDT_treat ///
		L_ttreat post co_code1 A_EDT B_LL B_female C_EDT_LL C_female
	
	gen A_EDT_elig = (A_EDT==1 & d_lleligible==1)
	gen A_EDT_no_elig = (A_EDT==1 & d_lleligible==0)
	gen nobt_ne_nolw = EDT_treat==0 & d_lleligible==0 & L_ttreat==0
	
	gen A_EDT_elig_fem=A_EDT_elig*female
	gen A_EDT_no_elig_fem=A_EDT_no_elig*female
	gen nobt_ne_nolw_fem = nobt_ne_nolw*female
	gen EDT_female=EDT_treat*female
	
	local yvar "inc_asset_aggregate3_fu co_outcomes_aggregate_fu outlook_life_aggregate_fu hh_decision_pca_fu"
	local xvars "female d_lleligible bl_bus_owner fu_dummy*"

*Panel A
	foreach var in inc_asset_aggregate3 outlook_life_aggregate hh_decision_pca {
		reg `var'_fu `var'_bl  EDT_treat female d_lleligible fu_dummy* bl_bus_owner, cluster(co_code1)
		gen xbhat = _b[_cons] + _b[`var'_bl]*`var'_bl +  _b[female]*female + ///
			_b[d_lleligible]*d_lleligible + _b[fu_dummy1]*fu_dummy1 + ///
			_b[fu_dummy2]*fu_dummy2 + _b[fu_dummy3]*fu_dummy3 + ///
			_b[fu_dummy5]*fu_dummy5 + _b[fu_dummy6]*fu_dummy6 + ///
			_b[bl_bus_owner ]*bl_bus_owner
		gen `var'_res = `var'_fu - xbhat
		reg `var'_res EDT_treat
		drop xbhat
	}

	reg  co_outcomes_aggregate_fu EDT_treat female d_lleligible fu_dummy* bl_bus_owner, cluster(co_code1)
	gen xbhat = _b[_cons] +  _b[female]*female + _b[d_lleligible]*d_lleligible + ///
		_b[fu_dummy1]*fu_dummy1 + _b[fu_dummy2]*fu_dummy2 + _b[fu_dummy3]*fu_dummy3 + ///
		_b[fu_dummy5]*fu_dummy5 + _b[fu_dummy6]*fu_dummy6 + _b[bl_bus_owner ]*bl_bus_owner
	gen co_outcomes_aggregate_res = co_outcomes_aggregate_fu - xbhat
	reg co_outcomes_aggregate_res  EDT_treat
	drop xbhat
	
	mhtexp inc_asset_aggregate3_res co_outcomes_aggregate_res outlook_life_aggregate_res hh_decision_pca_res, treatment(EDT_treat)
	drop inc_asset_aggregate3_res co_outcomes_aggregate_res outlook_life_aggregate_res hh_decision_pca_res

*Panel B
	foreach var in inc_asset_aggregate3 outlook_life_aggregate hh_decision_pca {
		reg  `var'_fu `var'_bl  EDT_treat EDT_female female d_lleligible fu_dummy* bl_bus_owner, cluster(co_code1)
		gen xbhat = _b[_cons] + _b[`var'_bl]*`var'_bl +  _b[female]*female + ///
			_b[d_lleligible]*d_lleligible + _b[fu_dummy1]*fu_dummy1 + ///
			_b[fu_dummy2]*fu_dummy2 + _b[fu_dummy3]*fu_dummy3 + ///
			_b[fu_dummy5]*fu_dummy5 + _b[fu_dummy6]*fu_dummy6 + ///
			_b[bl_bus_owner ]*bl_bus_owner
		gen `var'_res = `var'_fu - xbhat
		reg `var'_res EDT_treat EDT_female
		drop xbhat
	}

	reg  co_outcomes_aggregate_fu EDT_treat EDT_female female d_lleligible fu_dummy* bl_bus_owner, cluster(co_code1)
	gen xbhat = _b[_cons] +  _b[female]*female + _b[d_lleligible]*d_lleligible + ///
		_b[fu_dummy1]*fu_dummy1 + _b[fu_dummy2]*fu_dummy2 + _b[fu_dummy3]*fu_dummy3 + ///
		_b[fu_dummy5]*fu_dummy5 + _b[fu_dummy6]*fu_dummy6 + _b[bl_bus_owner ]*bl_bus_owner
	gen co_outcomes_aggregate_res = co_outcomes_aggregate_fu - xbhat
	reg co_outcomes_aggregate_res  EDT_treat EDT_female
	drop xbhat
	
	recode female 0=2
	mhtexp inc_asset_aggregate3_res co_outcomes_aggregate_res outlook_life_aggregate_res hh_decision_pca_res, treatment(EDT_treat) subgroup(female)
	drop inc_asset_aggregate3_res co_outcomes_aggregate_res outlook_life_aggregate_res hh_decision_pca_res

*Panel C
	foreach var in inc_asset_aggregate3 outlook_life_aggregate hh_decision_pca {
		reg  `var'_fu `var'_bl  A_EDT_elig A_EDT_no_elig B_LL C_EDT_LL nobt_ne_nolw female fu_dummy* bl_bus_owner, cluster(co_code1)
		gen xbhat = _b[_cons] + _b[`var'_bl]*`var'_bl +  _b[female]*female + ///
			_b[fu_dummy1]*fu_dummy1 + _b[fu_dummy2]*fu_dummy2 + ///
			_b[fu_dummy3]*fu_dummy3 + _b[fu_dummy5]*fu_dummy5 + ///
			_b[fu_dummy6]*fu_dummy6 + _b[bl_bus_owner ]*bl_bus_owner
		gen `var'_res = `var'_fu - xbhat
		reg `var'_res A_EDT_elig A_EDT_no_elig B_LL C_EDT_LL nobt_ne_nolw 
		drop xbhat
	}

	reg  co_outcomes_aggregate_fu A_EDT_elig A_EDT_no_elig B_LL C_EDT_LL nobt_ne_nolw female fu_dummy* bl_bus_owner, cluster(co_code1)
	gen xbhat = _b[_cons] +  _b[female]*female + _b[fu_dummy1]*fu_dummy1 + ///
		_b[fu_dummy2]*fu_dummy2 + _b[fu_dummy3]*fu_dummy3 + ///
		_b[fu_dummy5]*fu_dummy5 + _b[fu_dummy6]*fu_dummy6 + ///
		_b[bl_bus_owner ]*bl_bus_owner
	gen co_outcomes_aggregate_res = co_outcomes_aggregate_fu - xbhat
	reg co_outcomes_aggregate_res A_EDT_elig A_EDT_no_elig B_LL C_EDT_LL nobt_ne_nolw
	drop xbhat

	tab EDT_treat L_ttreat,m
	gen treat=.
	replace treat=1 if A_EDT==1 
	replace treat=2 if B_LL
	replace treat=3 if C_EDT_LL
	replace treat=0 if EDT_treat==0 & L_ttreat==0
	
	tab d_lleligible female,m
	
	gen group=.
	replace group=1 if female==0 & d_lleligible==1
	replace group=2 if female==0 & d_lleligible==0
	replace group=3 if female==1 & d_lleligible==1
	replace group=4 if female==1 & d_lleligible==0

	recode d_lleligible 0=2
	mhtexp inc_asset_aggregate3_res co_outcomes_aggregate_res outlook_life_aggregate_res hh_decision_pca_res, treatment(treat) subgroup(d_lleligible)
	drop inc_asset_aggregate3_res co_outcomes_aggregate_res outlook_life_aggregate_res hh_decision_pca_res

*Panel D
	foreach var in inc_asset_aggregate3 outlook_life_aggregate hh_decision_pca {
		reg  `var'_fu `var'_bl  A_EDT_elig  A_EDT_elig_fem A_EDT_no_elig ///
			A_EDT_no_elig_fem B_LL B_female C_EDT_LL C_female nobt_ne_nolw ///
			nobt_ne_nolw_fem female fu_dummy* bl_bus_owner, cluster(co_code1)
		gen xbhat = _b[_cons] + _b[`var'_bl]*`var'_bl +  _b[female]*female + ///
			_b[fu_dummy1]*fu_dummy1 + _b[fu_dummy2]*fu_dummy2 + ///
			_b[fu_dummy3]*fu_dummy3 + _b[fu_dummy5]*fu_dummy5 + ///
			_b[fu_dummy6]*fu_dummy6 + _b[bl_bus_owner ]*bl_bus_owner
		gen `var'_res = `var'_fu - xbhat
		reg `var'_res A_EDT_elig  A_EDT_elig_fem A_EDT_no_elig A_EDT_no_elig_fem ///
			B_LL B_female C_EDT_LL C_female nobt_ne_nolw nobt_ne_nolw_fem
		drop xbhat
	}
	
	reg  co_outcomes_aggregate_fu A_EDT_elig  A_EDT_elig_fem A_EDT_no_elig ///
		A_EDT_no_elig_fem B_LL B_female C_EDT_LL C_female nobt_ne_nolw nobt_ne_nolw_fem ///
		female fu_dummy* bl_bus_owner, cluster(co_code1)
	gen xbhat = _b[_cons] +  _b[female]*female + _b[fu_dummy1]*fu_dummy1 + ///
		_b[fu_dummy2]*fu_dummy2 + _b[fu_dummy3]*fu_dummy3 + ///
		_b[fu_dummy5]*fu_dummy5 + _b[fu_dummy6]*fu_dummy6 + ///
		_b[bl_bus_owner ]*bl_bus_owner
	gen co_outcomes_aggregate_res = co_outcomes_aggregate_fu - xbhat
	reg co_outcomes_aggregate_res A_EDT_elig  A_EDT_elig_fem A_EDT_no_elig ///
		A_EDT_no_elig_fem B_LL B_female C_EDT_LL C_female nobt_ne_nolw nobt_ne_nolw_fem
	
	mhtexp inc_asset_aggregate3_res co_outcomes_aggregate_res outlook_life_aggregate_res hh_decision_pca_res, treatment(treat) subgroup(group)
